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Abstract 

We present a systematic construction of the probability-current operator, based on a momentum 
power expansion of effective Hamiltonians. The result is valid in the presence of a Rashba term 
and when a D'yakonov-Perel contribution is included. We propose practical tools for spin-orbit 
engineering of semiconductor heterostructures. We apply this formalism to a paradigmatic system, 
the interface between two semi-infinite media, on one side a free-electron-like material and on 
the other side a barrier material with spin-orbit interaction. We show that the usual boundary 
conditions, namely the continuity of the envelope function and of a velocity at the interface, 
according to the BenDaniel-Duke approach, comply with the conservation of the probability current 
only when first- (Rashba-like) and second-order (free-electron-like) terms are taken into account 
in the Hamiltonian. We revisit the boundary conditions and we prove that the envelope function 
may be discontinuous at the interface. 

PACS numbers: 72.25.Dc, 71.70.Ej, 73.40.Gk 
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I. INTRODUCTION 



The probability current is a fundamental concept in quantum mechanics, which connects 
the wave-like description of a quasi-particle to the notion of transport current. When we 
consider a general Schrodinger problem where the Hamiltonian is 

fio = £ + ilW, (1) 

where the real potential il(r) is periodic in a crystalline solid and m is the free-electron 
mass, we are led to the usual definition of free-electron current probability: 1 



J f U} = Re 



m 



h 

— Im U* Vib] . (2) 
m 



However, in condensed-matter systems in the presence of Spin-Orbit Interaction (SOI), 
the potential is no longer real so that a redefinition of this quantity is mandatory. A 
debated example of this subtle point is provided by semiconductors-based systems, whose 
proper treatment requires consideration of the Hamiltonian 

H = H + Hso (3) 

with 

Hso = ^ (Vil x p) ■ a. (4) 

Following the arguments developed by Hoai Nguyen et al.,— it is reasonable to express 
the full Hamiltonian, involving SOI terms, as an effective Hamiltonian which consists of 
momentum-operator p-power series expansion: Indeed, beside the kinetic energy, quadratic 
in p, the SOI provides leading terms that are linear and cubic in p, known respectively as 
Rashba 3 and D'yakonov-Perel (DP) 4 terms. Then, since the SOI potential is non real, a more 
general definition of the probability current J[ip] has necessarily to be taken into account. 
Considering interactions that include higher-order polynomial terms in the Hamiltonian, we 
have to deal with an effective Hamiltonian of order n. 

Furthermore, an open question, strictly linked to the one above, concerns the definition 
of spin current (SC). Indeed, in semiconductor physics, that provides paradigmatic systems 
for spintronics, it is known that the SC standard definition, used by many authors,- - — can 
be suitably applied to two dimensional (2D) systems with Rashba SOI, but fails to describe 
spin-dependent transport phenomena in bulk cubic semiconductors, where SOI induces a 



DP term in the conduction band. The existence of extra-current terms was pointed out in 



Ref. 
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Drouhin et al— have shown that a redefinition of SC is mandatory to obtain a 
unified treatment, enlightening the fact that a properly-symmetrized spin-current operator 
J^m^], where t (I) refers to up (down) spin channel, gives unexpected results when applied 
to tunneling through evanescent states in GaAs barriers. 

As pointed out by Rashba in Ref. O, there are still concerns relying on the fact that a 
consistent theory of spin transport currents has not been formulated yet. From a general 
point of view, it means that we cannot immediately approach such a topic in terms of 
non equilibrium thermodynamics. In fact, a difficulty relies on the definition of system in 
order to formulate relevant balance equations and also on the boundary terms which should 
possibly be included in the effective Hamiltonian. Recently, Shi et alM have proposed an 
alternative spin-current operator, satisfying the continuity equation, that allegedly supports 
important conclusions concerning conservation of spin currents,—"— but which appears to 
rely on non-explicit assumptions (see Sec. HT1) . 

The inclusion of SOI in the Hamiltonian of a system has direct and practical consequences 
in heterostructures, where a consistent analysis of the tunneling phenomena is required. 
The pragmatic BenDaniel-Duke (BDD) approach, 17 that perfectly works when dealing with 
quadratic Hamiltonians under effective-mass approximation, cannot be straightforwardly 
extended because it is not always possible to ensure both the continuity of the envelope 
function and the conservation of the probability current, which is mandatory under steady- 
state conditions. Then it is necessary to revisit both the probability-current expression and 
the boundary conditions. This is in line with the ideas of Harrison.— 

In this paper, we present a systematic construction of the probability-current operator 
J, based on an effective Hamiltonian written as a p-power series expansion. We show 
the relation between the velocity operator and the current operator, evidencing the simple 
structure of the extra terms. This yields easy and compact calculations whereas explicit 
treatments in particular cases resulted in lengthy calculations.— The current operator can 
be subsequently used to build the SC operators according to the procedure described in Ref. 
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Then, we introduce proper matching conditions at the boundaries, which generalize 
the BDD procedure, the simplest efficient way to deal with semiconductor heterostructures. 
Finally, we illustrate our method on three examples: the case of a quadratic Hamiltonian, 
where we recover the usual situation (continuity of the envelope function and of the velocity), 
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the case where a Rashba term is added as a perturbation to the BDD Hamiltonian (there 
we find that the envelope function is continuous, but its derivative is discontinuous), and 
the case where a cubic DP term is added to the BDD Hamiltonian (where we prove that 
the envelope function cannot be continuous). 

The layout of this paper is as follows: Sec. [Ill we gi ye a general construction of current 
operators and a derivation of local properties. Sec. Ill, we introduce a general Hamiltonian 
as a n th -degree homogenous function of momentum-operator coordinates; we consis- 
tently derive the velocity operator and we show that a proper symmetrization yields the 
Hermitian current operator J. Sec. IIV[ we propose boundary conditions which are suitable 
to deal with heterostructures. Sec. [V] we apply our method to electron tunneling through 
a [110]-oriented GaAs barrier. Sec. IVIl we show how to extend the construction procedure 
to the spin current operators. 



II. GENERAL DEFINITION OF CURRENT OPERATORS 



A difficulty, that arises when a current operator is taken into account, relies on the correct 
definition of the system and of its boundaries: in fact, considering the density p of a physical 
quantity, we need to satisfy the continuity equation of J, defining a source term G, so that: 

| + V.J^ G . (5) 



As pointed out by Shi et al. in Ref. [13|, the continuity of J can be ensured by introducing 
a general source term G, as in Eq. El but the source term is not uniquely defined and 
this leads to possible confusion when considering the conservation laws in terms of non- 
equilibrium thermodynamic equations.— In any case, we first need to state clearly the local 
properties of a current operator, postponing the analysis of its global properties. For this 
purpose, we consider a linear operator A that does not explicitly depend on time and acts 
over a generic state ip. In the following we adopt the notation (^Cj = Aip^j = tp^A-ip 



used in Ref. |2|. The general Schrodinger problem reads: 

ih^ = Ui> (6) 

where "H may be any Hamiltonian. For example H may be equal to H (defined in Eq. [3]) or 
to H e fj (defined below in Eq. [26]). We explicitly develop the derivative of A with respect 
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to time: 

dt V7 dt 
and with the help of Eq. [6] we obtain: 



'A)=t(^Aij 



9 W)A1> + tfA(^ 



dt 



-— (Hil>) a^ + —^a(u^\ = 4 

in \ J in \ J %h 



(?) 



(8) 



If A is an Hermitian matrix (the elements of which are complex numbers, not differential 
operators) 

(% i>J f Aip = ^ AH * , (9) 

so that we can rewrite Eq. [S] in a more suitable way that we refer to as the local form of 
Ehrenfest theorem: 

a\{ A ) = \**{f%w). (io) 

The integration over the whole space leads to the well known Ehrenfest's theorem, whose 
global form is valid for any (possibly differential) Hermitian operator A: 



0_ 

dt 

We can write 



.4 



1 

ih . 



au 



tI>)-{W 



A 



i h 



am 



(ii) 



|(l)4lm(^{l,H}^) + ilm(^t 



An 



(12) 



with |a,6[ =ab + ba 1 and, by integration over the whole space, we get 



d 3 rlm (V = 0. 



(13) 



The time derivative of yAj is composed of two parts, concerning two different physical 
processes: we respectively recognize in Eq. [12] the divergence of the current and the source 
term G associated to A 



V- J, 



^ Im (V {1, V>) =-ilm(v t {A^-w}v) , (14) 



where any real potential U vanishes when taking the imaginary part of the anticommutator, 
and 



G = - Ira 

h 



A,n 



(15) 
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The above procedure has two advantages: first, we have expressed in a general form all the 
quantities entering Eq. [5] through commutators and anticommutators; then we have related 
the probability-current expression directly to the local properties of its corresponding oper- 
ator, without taking into account a closed system (such a procedure does not automatically 
imply that the integral of V • J a over the crystal only is zero). The choice of considering 
open systems makes the current operator involve Dirac distributions to deal properly with 
possible discontinuities at the boundaries of a subsystem. It has to be noted that it is always 
possible to include the source G term in the form of a current Jg, G = V • Jg so that the 
conservation equation becomes 

|(l) + V.(J,- Jo ) = |(i) + V.^0 (16) 

where J = J a — Jg is divergence-free in steady-state regime. For instance, if we look for 
Jg = V£/g, the potential Ug is a solution of the Laplacian problem AUg = G. Moreover, 
adding to Jg the term Vx^4 G , where Ag is an arbitrary vector field, does not affect the 
conservation equation. At this stage, the boundary conditions are not under control. Shi et 
al.— observe that it might often happen that 

I d 3 r G = (17) 
Jy 

where the integration is performed over the volume of the system (V). Then 

f d 3 r G = / d 3 rV-J G = / J G ■ ds = (18) 
Jy Jy Js 

where the volume integral is changed into a surface integral through Ostrogradski's theorem 

(here § is the surface limiting V and ds is the surface element along the normal to S). Such a 

relation is obviously satisfied provided that J^-ds = 0, i.e., provided that Jg is a tangential 

vector to §, which is "physically" reasonable. Shi et al. further assume that Jq "is a material 

property that should vanish outside the sample" : this is a more restrictive and questionable 

hypothesis. For instance in the case of a magnetic field, the effect of the associated vector 

potential cannot a priori be overlooked outside the sample. Anyway, let us assume that 

Jg = at the surface S. Following the calculation by Shi et al., it is straightforward to 

show, after partial integration where the boundary contribution cancels, that 

/ d ^K^ + ^ + ^r) = -/ dW - (19) 
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where Jg,x, Jc,y, and Jq, z are the Cartesian components of Jg- Then 



d 3 r J G = - / d 3 r r V ■ J G = - / d 3 r tG 



-— I d 3 r r Im(^ 



— I d 3 r Im (ip^r 



4> 



It is easy to check that, provided that 



A,r 



Ar,H 



ihvA, 



where 



ihv. Thus 



d 3 r 3 G = -- I d 3 r Im ( ^ 



Ar,H 



h 



d 3 r Im ( ift 



Ar,U 



A + ( d 3 r 3 A . 



Here, J a is the canonical current defined as 



v A -\- A v 
J A = Re (ijj^vAif)) = ^ ip. 



According to Eq. [HJ we can write 



d 3 r J 



G 



d 3 r 



d At 



dt 



d 3 r 



d I At 



dt 



Shi et al. define the effective current density as J G 

d (Av 



>G 



dt 



(20) 



(21) 



(22) 



(23) 



(24) 



We have the two following relations which define respectively the total current J and the 
effective total current J 



J=3 



A — JG, 



_ d (^4r 

J=^a — Jg — — 7 + ( Ja — J a 



(25a) 
(25b) 



Provided 3 a — J a = 0, i.e. when making the confusion between the canonical and the 
true currents (which is justified only for Hamiltonians up to second order in p, see Sec. 
Ill), the effective total current becomes J =d( Arj /dt, which is Eq. 5 in the papers by 



Shi et al.— and also by Zhang et al.,— and is the cornerstone of their further calculations. 
After a careful analysis, this relation appears to be derived under very special conditions 
so that it cannot be general. Moreover, the meaning of the so-called effective currents and 
their relationship with the true currents are not clear. Their use to tackle local transport 
equations is not justified. 



III. PROBABILITY CURRENT OF AN EFFECTIVE HAMILTONIAN 
A. Formulation of the general n th -order Hamiltonian 

Considering effective Hamiltonians, we deal with general expressions given by momentum 
series expansions, i.e., constructed from the energy expressed as wave- vector-component 
series expansion after the substitution {k — > — iV}. We write the effective Hamiltonian 
H e ff as follows: 

H eff = H P + V (r) (26) 

where V (r) is a potential which may be the potential of a single barrier or the one of a 
superlattice, for example, H p is such that 

H P = Y1 Yl c Ki)m,---An)PiW---hn) = J2H {n) (27) 

Tl l(k)£{x,y,z} n 

k=l,...,n 

where pi^) is the momentum operator associated to the l(k) Cartesian coordinate and where 
Q(i), ...,/( n ) are Hermitian matrices which are invariant under permutation of the subscripts. 
The abstract form of Eq. [27] allows us to perform easy calculations. In Sec. IIIIBI we show 
how to handle such a general expression to deal with concrete situations. 
Formally, we perform the identification 

c x ...c x Cy...c y c z ...c z = Cx...x,y...y,z...z 

(28) 

a p 7 " 7 

where a, (3, and 7 are integers. We obtain 

H {n) = (c x p x + c y p y + c z pzf ■ (29) 

Given Eqs. |27J [281 an d EU let us note that only terms such as c xx or c xy (for n = 2) are 
meaningful, a term such as c x being only a trick in the calculation. 
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Alternatively, one can write 



H {n) = caPl P^PyP2 (3°) 

a+/3+7=ra 



with 



r <*/3 7 _ n - a J3 7 foi \ 

° ~ a!/3! 7 ! C ^ ^ 
We are now in a position to tackle the problem of velocity, first when the Hamiltonian 
H takes into account the SOI, and, second, when the Hamiltonian H e ff is an effective 
Hamiltonian. 



B. Velocity operator in presence of SOI interaction 



It is usually admitted that the velocity operator v is equal to &H/dp whatever the 
Hamiltonian H. However, to the best of our knowledge, the derivation can be found only 
when H is quadratic in p. Therefore a general derivation, in particular in the case of 
effective Hamiltonians, is mandatory We start from Ehrenfest's theorem (valid whatever 
the Hamiltonian %) 



(v) 



djr) 
dt 



H, r 



(32) 



if (i/h) 



n, ? 



dH /dp), then v = dH / dp because two linear operators which have 



the same mean values over all possible states are equal: ( A 



B) implies that A — B. 



20 



Practically, it is enough to show that {i/h) 



H, r 



&H/dp to prove that v = &H/dp. 

Ho, r 



First, considering the case "H = H which contains a term (Eq. [T]), (i/h) 
(h/im) p= dHo/dp, for a system described by an Hamiltonian quadratic versus momentum 
components, and we obtain the velocity v : 

- dH 



v 



dp 



(33) 



Second, we have to check that this relation still holds in the presence of SOI where the 
Hamiltonian is H = H = H + Hso (Eq. [3]). In other words we want to show that 



dH 
dp ' 



(34) 



We know that vq = [if h) 



H , r 



dHo/dp. To show that Eq. [34]is valid, it is enough to 



show that vso = dH so /dp, which will give v= dH /dp with v= vq+v^o- A straightforward 



calculation yields 

which proves Eq. [34j the derivative of the Hamiltonian, with respect to the momentum 
operator, still provides a suitable definition of the velocity when the SOI term is taken into 
account. 



C. Velocity operator with an effective Hamiltonian H e tt 



We generalize the results obtained in Sec. IIIIBl to the case of a generic effective Hamil- 
tonian H = H e ff. Again we exploit Ehrenfest's theorem, as written in Eq. [321 Considering 
for instance the x component, we verify that 

hdH^ 



H {n) 



X 



or 



i 

h 



a+/3+7=n 



H p , r 



% j % 



h 



which proves that 



dH e j f 
dp 



OH, 



eff 



dp 



(36) 



(37) 



(38) 



Using Eqs. |28]|3~T| it is then easy to calculate the j component (j = x, y, z) of the 
velocity operator associated to H^: 



■An) 



QH(n) 



ncj (c x p x + CyPy + c z pz 



\n— 1 



(39) 



We introduce the scalar product between the momentum p and the velocity operator v 



(n) 



%v x n) + P y v y n) + p# = n [c x p x + c y p y + c z p z ) n = nH^. 



(40) 



With this notation, vo, introduced in the paragraph IIII Bl is such that vo = v^ 2 ^. Eq. HQ] 
means that 

p • v (n) = nH {n) (41) 

and eventually 

" ' ^ (12) 



H ef ^ = (p ■ ^ ^ (n) j il> + V1> = Eip. 
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As pointed out in Sec. [TTJ we are allowed to define current operators in open systems 
provided that we properly take into account their boundary conditions. In Appendix [Bj we 
show that performing the proper symmetrization according to the following rule (Eq. |4"4"|) . 
we find a probability current operator that for the j^-Cartesian component reads: 

Z (ro) = £ J i n) ( r o) (43) 

n 

Jj (r ) = £ C 3,m,-Mn-1) [ S r Pl(l)-Pl(n-l) + Pl(l)$T -Pl(n-l) + ••••+ Pl(l)-Pl(n-l)$T ] 

l(k)£{x,y,z} 

fc=l,..,n— 1 

(44) 

where S ro = 5(r — ro) is the Dirac distribution. We must still verify that the divergence 
of the current, calculated with the operator defined by Eq. HU satisfies the conservation 
equation for the density of probability (Eq. [T4lwhen A is the identity). It is straightforward 
to show (see Appendix |B]) that the divergence of the probability current can be written as: 
2 

VJ M = £ V ' jW M = Im E £ Yl C m),-,l(n-l) (V> |PiPJ(l)-P/(n-l)| V>) 

n n J = {x,1/,2} !(fc)g{a:, ! ;,z} 

fe=l,..,n-l 

n (45) 

where we again use the notations Ai[)j = tp^Aip as in Ref. M. Then, we recover all 
the terms of Eq. [TH so that Eq. HH indeed provides a correct definition of the current 
operator. Obviously, adding a term proportional to the curl of any vector field would not 
affect the result. Such a definition of J provides an unambiguous and general tool for 
evaluating the probability current. Provided the Hamiltonian of the whole system is known, 
this probability-current operator guarantees the requirements of the continuity equation. 
Now it is useful to introduce the Hermitian symmetrized velocity operator 

^"Vo) = ^E C m),-,Kn-l) [KPl{l)-Pl{n-l) +Pl(l)---Pl(n-l)8r ] (46) 
fe=l,..,n-l 

For example for n > 2, the comparison between Eqs. HH and H6] clearly shows 
that jj n \r ) contains n — 2 extra terms, which are straightforwardly obtained from 
dH e ff/dp. For instance, with H e ff = p", we have dH e ff/dp = np n ~ 1 , so that t)( n )(r ) = 
(n/2)(5 ro p n - 1 + p n - 1 5 ro ), whereas J (n) (r ) = (c^oP"" 1 + pKp"' 2 + ••• + F"" 1 ^)- As 
shown in Ref. extra terms are specially important for evanescent waves. Therefore, in the 
following we deal with tunneling problems. 
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IV. BenDaniel-Duke-LIKE FORMULATION AND BOUNDARY CONDITIONS 



We stress that the central question when defining the current operators and related quan- 
tities is the proper definition of the system and of its boundaries. Dealing with heterojunc- 
tions, where each bulk medium is described by the relevant Hamiltonian, requires defining 
proper matching conditions at the boundaries. In this sense, the BDD Hamiltonian^ is the 
simplest smart approach that allows solution of the Schrodinger equation over the whole 
space while it guaranties the conservation of the probability current at the interface. The 
principle is the following. Let us consider a one-dimensional problem and two different media 
for x < and x > 0. Each medium is characterized by its own Hamiltonian. The question 
is to find a solution of the Schrodinger equation, made of eigenfunctions of the relevant 
band of the two bulk materials, which ensures the continuity of the probability current at 
the origin. In this sense, the problem is analogous to a scattering problem, where the wave 
functions are determined only at some distance of the scattering potential. Proper matching 
conditions relevant to the extension of the bulk envelope functions at the origin will allow 
one to determine the envelope function over the whole space. For that, BBD propose writing 
an Hamiltonian over the whole space as p x [1 /2m (x)} p x + V (x) where m (x) is the effective 
mass in each medium. The integration of this BDD Hamiltonian around the boundary au- 
tomatically ensures the continuity of the probability current of Eq. provided that ip (x) 
and [l/m(x)] [dip/dx] are continuous. 

Now, consider two regions (1) and (2) and assume that each region is made of a given 
crystalline material. We look for the envelope function, solution of the Schrodinger equation, 
which is made from plane waves which are eigenstates of the crystal, inside each material. 
Observe that, near the interface, the crystal periodicity is broken so that the true Hamil- 
tonian and the true eigenfunctions will become involved. The principle is then to define 
proper matching conditions applying to the prolongation of the envelope function at the 
origin. For that purpose, we consider a volume V, limited by a surface §, that surrounds an 
interface portion. Similarly to the BDD technique, we start from Eq. H2] and we integrate 
the Schrodinger equation over V. Using Ostrogradski's theorem, when V tends to zero, we 
obtain 




(47) 



where ds is normal to the surface §. 
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For a one dimensional case with the interface at the origin, Eq. 1471 becomes: 



lim 



y i«<»ty 

^ n 



+e 



(48) 



Let us again emphasize that no information is obtained on the true wave function near 
the origin. Eq. [48] does not ensure either the continuity of the envelope function or the 
existence of derivatives at the interface. 

As an illustration, let us consider the case of a Rashba Hamiltonian 



H eff = ap + bp 2 (49) 

where a and b are two Hermitian matrices. According to Eq. HH1 we can write down the 
first continuity condition as follows: 

[aip + bpip}t £ £ = 0. (50) 

Using this condition to solve the problem, and adding a priori the continuity of the 
envelope function at the interface as a second condition, we verify that the probability 
current is indeed continuous at the interface: 

J [^] = (tf | a + bp | ip) + c.c. (51) 
Then, the jump of the derivative of the wave function at the interface is determined by 

[atim- (52) 



It is clear then that the BDD approach, introduced to solve a problem with a quadratic 
Hamiltonian, is also suitable to obtain a solution when a Rashba contribution is added; then 
we can say that up to the second order in the momentum-power series expansion of the 
Hamiltonian, the continuity of a "generalized velocity" (see Eq. l4~8l) and the continuity of 
the wave function at the interface imply the conservation of the probability current at this 
point. Remarkably, the boundary conditions that we need to solve the problem drastically 
change when moving to the case of a DP Hamiltonian with cubic terms. The crucial point, 
that we address in the following, is that we cannot make any hypothesis about the continuity 
of the wave function because, if we need to ensure probability-current conservation at an 
interface, we must accept an envelope function ip which is no longer continuous. 
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To give an insight into the expression of the current operator and into the conservation 
of the probability current, let us again come back to an interface between two semi-infinite 
one-dimensional media (1) and (2). In each bulk crystal, the relevant Hamiltonian is 

H r = J2^ n) + V r (53) 

n 

with 

= 7 Wf (54) 

with r = 1 or r = 2 depending on wether x < or x > 0. H r admits the eigenf unctions (p r , 
associated to the fixed energy E which verify 

H r ip r = Etp r . (55) 

Near the heterojunction, the spatial periodicity is broken, so that over a few Wigner-Seitz 
cells, the electron states are no longer pure Bloch states. We consider two coordinates, —w\ 
and w 2 , so that, in the regions ]— oo, — Wi] and [w 2 , +oo[ the electronic structure remains 
unaffected. In these regions, the Hamiltonian writes 

$ = Q(- x - Wl )H 1 + Q(x-w 2 )H 2 (56) 

where 6 (x) is the Heaviside function. In the vicinity of the heterojunction, ]— wi,w 2 [, 
the form of the Hamiltonian and of the wave functions are not explicitly known. We con- 
sider a wave function ip which is an eigenstate of the Hamiltonian over the whole space 
at energy E. We expect that, over the domain ]— oo, — wi] U [u>2,+oo[, ip coincides with 
\& = 6 (—x — Wi) ipi + 6 (x + w 2 ) <p 2 - We thus expect the following relation to be satisfied 



^6 (-x) + ip 2 & (x) f)V)=E [(</?! | v^i) (1) + (v^2 I V^2> (2) J (57) 

where (|)^ means summation over the bulk part of region (r). Observe that: 

H^m = ^pTQ (-x - wi) (f! + V X Q (-x - Wl ) <f! 
= Q{-x-w l )(^ i C ) r + V l ) Vl 

+ ih-y[ n) [5 (x + w^pT- 1 +p5{x + Wl )f l ~ 2 + ... + f l ~ l 5 (x + Wl j\ n 

= G(-x- w ± ) (i[ n) F + V^j <p x + iftjfVi = 6 (-x - w x ) E(p! + iftjfVi (58) 
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and similarly 



H^m = Q(x + w 2 ) E V2 - ih4 n) <p 2 



(59) 



So that 



(pi® (—x) + (^ 2 (x) 



(fi I Vi)(i) + (^2 I Va)(2) 
+ z/i(j* n) M]- J? n) [</?i(-wi)] 



This implies that 



(60) 
(61) 



jW [cp 2 (w 2 )} = J[ n) [cpt(- Wl )}. 

The important point is not the conservation relation, which might appear as physically 
obvious, but that, in Eqs. and E2U the symmetrized current operator is automatically 
generated in the form derived in Eq. H3J providing a physical insight into this mathematical 
expression. Because we only deal with wave functions at some distance from the heterojunc- 
tion, the continuity of the true wave function at x = does not implies the continuity of the 
envelope function \1/ which may be discontinuous. This is in line with the considerations of 
Harrison. 18 Eq. HH1 and ED] generate a set of boundary conditions relevant to the tunneling 
problem. 



V. THE [110]-ORIENTED GaAs BARRIER 

We analyze the case of electron tunneling under normal incidence through a GaAs [110]- 
oriented barrier, which was shown to be non trivial and solved in special cases in Ref. U- 
Hereafter, we apply the tools and boundary conditions presented in this paper to solve it in 
a more simple and general manner. We confirm and generalize the results derived in Ref. Q 
In particular, we are able to solve the problem of an heterojunction between a free-electron 
media and a semiconductor without inversion center, where the DP field is a step function, 
which remained puzzling. In the [110] direction, the DP Hamiltonian is: 

where + (— ) refers to the up (down)- spin channel quantized along the DP field direction. 
We consider as solution a general wavefunction written as follows: 

ij> = aUo+ jr^JcP^ e %xz (63) 
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where ipo is the zeroth order function that is a solution of the tunneling problem with energy 
E and with the potential V when SOI is turned off. Here 7 c /c 2 — E—V, a and /3 are complex 
parameters to be determined, and x is a rea l ( see below) wavevector component which is 
added to k when SOI is turned on. We have the relations 



pip = a \ pip + i 



3 



p 2 ip 



h 2 



E-V 



7c 



- x 2 ) 4> + iftxv^- 



We calculate the velocity operators from Eq. 



7c- , 7 -2 



2 3 

and, according to Eq. HH1 we find the matching condition 



(64) 
(65) 

(66) 



7c 



1 ± — X ) 

7c 



( 



E-V 



\ 7c 



-r up 



+e 



(67) 



which is a generalization of Eq. (3.50) of Ref. 

Now, we have to satisfy the conservation of the probability current 

E-V 



7c 

h 2 



1 ± — X ) + i> (pip)*] ± 



7c 



7c 



- X 



\ij\ 2 ±^L\p^ 



(68) 



We obtain here an important result: The envelope function cannot be continuous at the 
interface. Indeed, assume ip to be continuous. Then, after Eq. [671 we see that the last 
term in Eq. EHJ that we rewrite as ±(l/2h 3 ) (7/7^) 7c If0| 2 must be continuous. This is 
not possible since ■y c p4>o is continuous (unless 7/7^ is almost continuous, which would be 
fortuitous). 

We have to determine ip complying the boundary conditions, which is not simple because 
the expression providing the current is not a linear function of ip. However, if we consider 7 
as a first-order quantity and look for a solution to first order only, the result is surprisingly 
simple, as shown below. From the Schrodinger equation - Eq. [62]-, we find that x verifies 



7c {2k X + X 2 )±l {k 3 + 3k 2 X + 3k X 2 + X 3 ) = 



then 



X 



4 Tc 



=F 



1 7 E-V 

4 7c 7c 



(69) 



(70) 
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As stated above, \ is a real quantity. For each spin, there are two others roots of the 
cubic equation (Eq. |69|) which are much larger than the width of the Brillouin zone ; These 
two roots are of the order of 7 c /7 which is about 2 A (two times the Brillouin zone width) in 
GaAs (see Fig. 4 of Ref. [j) and have no physical meaning. Note that, the cubic DP term, 
obtained from perturbation expansion, only holds for small wave vectors, a few percent of 
the Brillouin zone, so that taking into account these two other roots would be meaningless. 
From Eq. [63] we see that, upon tunneling, the up- and down- spin electrons undergo opposite 
phase shifts, which is equivalent to a precession around the DP-field direction. This would 
be quite intuitive if the field were not a complex quantity, and constitutes a prediction which 
can be experimentally tested. Let us calculate the current at the interface J[ip (0)] to first 
order 



J[i/j (0)] = H 2 ^ [if)*P^o + (p Vb)1 + |: 



+ 



ry E - V 

2|a| 2 ( X -Im/3)± 



7c 

h 3 



7c 7c 

i i2 27 c 7 

— a 1m p =F 

1 1 E-V ^ + 2 7c 



IV>o| 



(71) 



where the values of ipo and of its derivative are calculated at the origin. Observe that with 
the choice 



|a| = 1 and Im/3 = — x 
the second and the last terms of Eq. [7U vanish so that 

J[*l> (0)] = |f bPoP^o + (P^o)*] = J f VI> (0)] 



(72) 



(73) 



where J^[ip (0)] results from the application of the free-electron current operator on ip (0). 
Thus, we obtain another essential result: To first order, turning on the SOI does not alter 
the value of the probability current. Consequently, to solve the problem we have only to 
show that ip, given by Eq. [63] and with the conditions defined in Eq. [721 can match the 
boundary condition expressed by Eq. (67J We obtain 

7 C [cap + H (ia{3 — x)} V'o — 7c [<xp + h [ot (i Re j3 + x) ~ x]\ i'o continuous. (74) 



The continuity of Eq. [74] can always be ensured by taking a = 1 and Re 
we calculate if) according to Eq. [63] and deduce the jump of if> at the interface 



X 



_h{E-V)_ 



o+ 



(7cP>o) = 



1_ 

lc J 0- 



0+ 



(7cP>o) • 



0. Then, 



(75) 
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VI. SPIN CURRENT 



We have shown in Sec. |V] that Eq. [44] provides a general and symmetrized definition of 



the probability-current operator. Following the conceptual scheme developed in Ref. |l2 



we 



can define the spin currents in the up- and down-spin channels by taking A = tt s , where n s 
is the orthogonal projector on the spin basis (s = ±). Then the SC current SJ u j [ip], that 
arises from the difference between the up-spin and the down-spin current, is obtained by 
taking A = <r u , the Pauli operator along the u direction defining the quantization axis. It 



is straightforward to see that, as in Ref. 



12 



, the j-component of the spin-current operator 



is obtained from the j-component of the probability-current operator after the substitution 

- 1 

-j,l(l),...,l(n) - 2 



C 'j,l(l),-An) = o{°"u, C,^!),...,^)}. (76) 



VII. CONCLUSION 



We have proposed a systematic procedure to construct properly-symmetrized current op- 
erators. We have extended the BenDaniel Duke approach™ to deal with heterostructures 
where SOI is included, introducing generalized boundary conditions, which allow us to con- 
sider open systems. We have shown that up to second order, usual matching conditions and 
in particular the continuity of the envelope function at an interface, yield solutions which 
comply with the conservation of the probability current. This no longer holds as soon as 
a cubic term is included (D'yakonov-PereF term). We have illustrated our findings on the 
model case of a [110]-oriented GaAs barrier, which had already been addressed in Ref. [2!. We 
recover and generalize the results of Ref. in a more simple and direct treatment. The tools 
we have developed can be applied to the holes in the valence bands or to the electrons in 
the conduction band so that they should be important for semiconductor-based spintronics. 
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Appendix A: Symmetry properties of current operators 



In Sec. |Vl Eq. [TOJ we derived the local form of the Ehrenfest theorem for a general 
operator A and deduced the expression of the associated current J A . First, consider the 
case where A = I, where / is the identity and the quadratic Hamiltonian p 2 /2m. We 
rewrite Eq. [TU1 as 

-V • J 



^H 2 = -V-Re(VP^ 
at \ m 

We recover the usual expression for the free-electron probability current 

J ty] = Re f ^ — tp 



m 



(Al) 



(A2) 



Note that: 



^ht) ~ GO 



i 

?7i 



(A3) 



where Kq is the time-reversal Kramers operator for a spinless particle, which consists of 
taking the complex conjugate in the r-representation. Let us check the expression of the 
current operators we defined under time inversion symmetry. For this purpose we consider 
the term 

-2ihV -3 A = 2ilm (V {l, #} ip) = ip ] AHip - (V AHi^* + ip ] H Aip - (V H AipJ 

(A4) 

First, look at the term ip^ A Hip 
(kip AHKip^j = (k ip WAHKip\ = (k ip W AK Hip^J = -s A (k ip 



e A K ip 



K AHip 



-e A [ip 



AHip 



R^ K AHip 
(A5) 



Here, K = RK is the Kramers operator for a particle with spin 1/2, R 
= R' 1 ^, and e A = ±1 depending wether A verifies^ 

KAK = e A A or R ] AR = e A A*. 

Similarly, for the term ip^ 1 H Aip 



-io„ 



(A6) 



Kip 



H AKip ) =-e A [K Q iP 
= -e A (k ip 



R? K H Aip 
k HAip) = -e A [iP 



HAip 



(A7) 
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Thus, we obtain 



2iIm(v t {A#}V>) = V» f \A,H]^ + e A (Ki?j {£h} (kA . (A8) 



We conclude that the general expression for the current of A is 

1 



(A9) 



Appendix B: Complete derivation of the current operator J 



We are interested in finding the form of the current operator J = (^J x , J y , J z ^j for 
an Hamiltonian + V (r) - the current operator being j( n ) - and more generally for 
the Hamiltonian H eff = H p + V (r) = £ n #W + V (r) (Eqs. I2E H2T|) - the current op- 
erator being J. For an Hamiltonian p 2 /2m + V (r), it is known 2 - 2 that the j th compo- 
nent of the current operator (j = x, y, or z) at the point r is of the shape jjP (r ) = 
(l/2m) [S ro Pj +pj 8 ro \; With the notation of Eqs. [2611271 = J2mz{*,y,*} Q(i),«(2)P/(i)Pi(2), 

fc=l,2 

jf ] ( r o) = J2i(i)={ x , y ,z} c m) [KPm +PiwK], Q(i),i(2) = (l/2m)fy iy(2 ). The aim of this 
appendix is to show that, for an Hamiltonian H^ n > + V (r), the following form of the j th 
component of the probability current operator 

Jj\r ) = C JAl),-Mn~l) [$r Pl(l)Pl{2)---Pl(n-l) 

l(k)e{x,y,z} 

fc=l,..,n-l 

+Pl(l)5r Q Pl(2)---Pl(n-l) + ... + P«(l)^(2)---K(n-l) 5 ro] ( B1 ) 

gives back Eq. [TH The Dirac distribution interacts with the mixed powers of the current 
operator so that the symmetrization procedure used in the construction of jj n ^(r ) provides 
(n — 2) further summations with respect to (r ). The two definitions coincide only up to 
n = 2. The extra terms are crucial in order to satisfy the continuity equation. We evaluate 
every term over a generic state for example the second term is of the shape 



(lf> \pl(l)$r Pl(2)---Pl(n-l)\ 1p) = I d 3 r 1p*pi { i)d ro Pi( 2 )-Pl(n-l)1p 

d 3 r (pi(i)ipy <y ro pi( 2 )...Pi(n-i)V' 



[Pi(i)tp (r )] Pi{2)-Pl{n-i)^ (ro) • (B2) 
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Then the j Cartesian component of probability current for ci generic stcite Jj \ib] can be 
written as: 



J 



(n) 



■0 



jj n) (ro) W) = ^2 C 3,l(l),-An) [^Pl(l)-Pl(n-1)^ + - 



l(k)e{x,y,z} 

fe=l,..,n-l 



+ {Pl{l)-Pl(k-l)^) Pl(k)-Pl(n-l)i> + - + {Pl(l)-Pl(n-l)ll>) i> (B3) 

where ip = ip (r ). From Eq. IB3[ we can find the generic divergence term related to the 
derivative with respect to pf 



(n) 



E 



Cj,/(1),...J(„) 



l(k)e{x,y,z} 

k=l,..n— 1 



+ (P/(l)-%-l)^) 7 PjPl(k)--Pl(n-l)^ ~ (pjPi(i)...pi( k -l)^) Pl(k)-Pl(n-1)^ 
+ {Pl(l)---Pl(k)^y PjPl(k+l)--Pl(n-l)^ ~ {PjPl(l)-Pl(k)^) Pl(k+l)-Pl(n-l)i> 

+ ... + {p l{l) ...pi {n ^ 1) i)) ] Pji) - (^(i)...^^.!)^) 1 ^ • (B4) 

In Eq. IB4I all the terms that have the same order in k (two consecutive terms but the 
first one and the last one) vanish after summation over j: 



E E <* 

J={x,J/,2} l(k)e{x,y,z} 

fc=l,..n-l 



+ (PKi)---^(fe)^) t P^(fc+i)--Pz(ri-i)V' =0 (B5) 
Then the only terms still remaining in the summation are: 



J2 ^ 



(n) 



P 



.J(n) 



j={x,y,z} 



Cj,Z(l),...,Q (n ) 



V' t PiP/(l)---Pi(™-l)^ ~ (PjPl(l)-Pl{n-l)^) i> 



E E 

j'={x,t;,2} l(k)£{x,y,z} 
fc=l,..,n-l 

j = {x,y,z} l(k)£{x,y,z} 

fc=l,..,n-l 



(B6) 



Now V • J (n) = (i/ft) P • J (n) [iff] and Eq. |B7] results in a collection of pure imaginary 
terms and the final expression for the divergence of the probability current reads: 

V • J (n) = -- Im ^ Y C j,i(i),-M) (V 1 \PjPl{i)-Pi{n-i)\^) ■ (B7) 



j = {x,y,z} l(k)S{x,y,z} 
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